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Abstract 

We analyse conditions for an evolution equation with a drift and fractional diffusion to 
have a Holder continuous solution. In case the diffusion is of order one or more, we obtain 
Holder estimates for the solution for any bounded drift. In the case when the diffusion is of 
order less than one, we require the drift to be a Holder continuous vector field in order to 
obtain the same type of regularity result. 

1 Introduction 

This papers concerns the solutions of an equation of the form 

u t + b ■ Vit + (-A) s ii = /. (1.1) 

We will prove that the solution u becomes immediately Holder continuous. In the case s £ 
[1/2, 1), all we require for the vector field b is that it is bounded. In the case s £ (0, 1/2), we 
require b to be in the Holder class C 1_2s . 

The case s = 1/2 is often referred to as the critical case because the diffusion is an operator 
of the same order as the advection. Both terms in (1.1) have the same weight at all scales and 
no perturbation techniques can be applied to obtain regularity of the solution u. This critical 
case was already studied in [11], where the result was used to establish the existence of classical 
solutions for the Hamilton- Jacobi equation with critical fractional diffusion. The main ideas in 
this paper originated in that work. In future work, we plan to apply the result of this paper to 
other nonlinear equations like for example conservations laws with fractional diffusion. 

It is interesting to compare this result in the critical case s = 1/2 with the result in [4]. In 
that paper the same equation is studied and the Holder continuity result is obtained using De 
Giorgi's approach for parabolic equations. The only difference is that the assumptions for b in [4] 
are that it belongs to the BMO class and it is divergence free. In this paper (or in [11]) we do not 
assume b to be divergence free, but we need the slightly stronger assumption that b £ L°° instead 
of just belonging to the BMO class. The methods used in the proofs are of a different nature. In 
[4], thanks to the assumption that divfr = 0, they use variational techniques (in the style of Dc 
Giorgi), whereas in this article our techniques are purely non- variational and are ultimately based 
on quantitative comparison principles which are arguably simpler. An alternative approach is used 
in [8] where the Holder continuity is obtained through the study of the dual equation. They still 
require that div b = in order for the dual equation to have the same advection-diffusion form. 

The case s < 1/2 is known as the supercitical case because the diffusion is of lower order than 
the advection. This means that in small scales the drift term will be stronger than the diffusion 
which seems to suggest that it would not be possible to prove any regularity for u. However, if we 
assume a Holder modulus of continuity for b, we can compensate this deterioration of our control 
of the drift term in small scales by a change of variables following the flow of b. The supercitical 
case is the main focus of this article. 
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It is interesting to compare our result in the case s < 1/2 with the result of [7]. In that article, 
Constantin and Wu explore the applicability of Caffarelli-Vasscur [4] techniques to the supercritical 
regime. They need to assume that b £ C 1_2s and divfe = 0. Therefore, in the supercitical regime, 
our non-variational method provides a strictly better result than the De Giorgi's approach, since 
the regularity assumption on b is the same, but we do not need to assume div 6 = 0. This result 
is used in [6] to obtain the well posedness of a modified supercitical quasi-geostrophic equation. 

The case s > 1/2 is often referred to as the subcritical case because the diffusion is of higher 
order than the advection in (1.1). The method in this paper provides a Holder continuity result 
for it if 6 £ L°° in the subcritical case. This condition is far from sharp, since one could obtain 
the same results under more general conditions. In [2], heat kernel estimates are obtained using 
perturbation techniques for s £ (1/2,1) for vector fields which belong just to a Kato class. In 
[12], Holder estimates are obtained in the case s £ (1/2, 1) but for b ■ V« replaced by a term with 
supcrlincar growth respect to Vu. The subcritical case is not the main objective of this paper. 
We still include it in our main theorem to show what result the method provides. 

Here is our main theorem. 

Theorem 1.1. Let u be a bounded function inW 1 x [0,1] which solves the equation (1.1). Assume 
that s £ (0,1), / £ L°°. Ifs £ [1/2,1), we assume b £ L°° . If s £ (0, 1/2), we assume b £ C 1 " 28 . 
Then u is Holder continuous for positive time t £ (0, 1]. Moreover, there is an estimate 



,\x~y\ a + \t-r\ a / ( - 2s '> 
t a /( 2s 



\u(x,t) - u(y,r)\ < yi Z (NIl- + 



for every x, y £ -B1/2 and < r < t < 1, where the constants C and a depend on s, n and ||6|| c i_2 S 
but not on u or f . 

It is our intention to keep this article as simple and accessible as possible. The fractional 
Laplacian could be replaced by arbitrary integro-diffcrential operators of order 2s, even with 
discontinuous coefficients, as in [11], but we chose to keep only the fractional Laplacian in order 
to make the article accessible to a larger audience. Moreover, we will concentrate on proving the 
a priori estimate for Theorem 1.1, assuming that we have a classical solution. In the last section, 
we will discuss the notions of weak solution and how the proof would adapt to those scenarios. 



2 Scaling 

If u satisfies the equation 

u t + b- Vu + (-A) s u = / 
then the function u x = A~ a u(\x, X 2s t) satisfies 

+ \ 2s ~ 1 b • Vu A + (-A) s m a = A 2s " Q / A 

with f x (x,t) = f(Xx,X 2s t). 

We define the parabolic cylinders in terms of the scaling of the equation: 

Q r :=B r x [-r 2s ,0]. 



3 The diminish of oscillation lemma 

The proof of the diminish of oscillation lemma uses essentially the same ideas as the proof in 
[11] for the case s = 1/2. The generalization to the full range s £ (0, 1) does not present major 
difficulties. In this section we present the proof in full detail. 
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Lemma 3.1 (point estimate). Let u < 1 in K™ x [—2,0] and assume it satisfies the following 
inequality in Br x [—2,0]. 

ut-A\Vu\ + (-A) s u<e (3.1) 
Where R = 2 + A. Assume also that 

\{u < 0} n (Si x [-2, -1])| > fi > 0. 

Then, if So is small enough there is a 9 > such that u < 1 — 9 in B\ x [—1, 0]. 
(ifte maximal value of Eq as well as the value of 9 depend only on s, fi and n) 

Remark 3.2. Note that if for some bounded vector field b, the function u satisfies the equation 

u t + b- Vu + (-A) s u = f, (3.2) 
then u also satisfies the inequality (3.1) with A = \\b\\ L<x> . 

Proof. Let m : [—2, 0] — > K be the function such that: 



m(-2) = 0, 
m'{t) = c Q \{x e B x : u(x,t) < 0}| - Cim{t). 

The solution of this ODE for m can be computed explicitly 

m(t)= [ c Q \{x :u(x,s) <0}nB 1 \ e - Cl{t - s) dt 



(3.3) 



We will show that if Co is small and C\ is large, then u < 1 — m(t) + 2eo in B\ x [—1,0]. Since 
for t e [-1,0], 

m(t) > c e- 2Cl \{u < 0} nBi x [-2,-l]| > c e~ 2Cl n. 

We set 9 — coe~ 2Cl n/2 for eq small and finish the proof of the Lemma. 

Let /3 : M — > R be a fixed smooth non increasing function such that /3(x) = 1 if x < 1 and 
/3(x) = if x > 2. 

Let r](x, t) = (3(\x\ + At) = (3(\x\ — A\t\). As a function of x, n(x, t) looks like a bump function 
for every fixed t. The bump shrinks at speed A. Recall that A will be the L°° bound for the 
vector field in the drift term, so n is a bump that is shrinking at the maximum possible speed. 

The strategy of the proof is to show that the function u stays below 1 — m{t)r](x, t) + e (2 + 1). 

Note that if there was no diffusion, 1 — mr](x,t) + e (2 + t) would be a supcrsolution of the 
equation for any fixed constant m. We would be able to prove that u < 1 — 9 in Qi only in the 
case that it(x, —2) < 1 — mf/(i, —2). The diffusion is the reason why we can use a variable value 
of m and obtain an actual improvement of oscillation. 

At those points where r\ = (precisely where \x\ > 2 — At = 2 + A\t\), (—A) s ri < 0. Since n 
is smooth, (— A) s ?7 is continuous and it remains negative for r\ small enough. Thus, there is some 
constant /3i such that (— A) s ry < where r\ < f}\. 

In order to arrive to a contradiction, we assume that rj(x,t) > 1 — m(t) + Eq(2 + t) for some 
point (x,t) G B\ x [—1,0]. We look at the maximum of the function 

w(x, t) = u(x, t) + m(t)ri(x, t) - e (2 + 1). 

Since we are assuming that there is one point in B\ x [—1, 0] where w(x, t) > 1, w must be larger 
than 1 at the point that realizes the maximum of w. Let (xq, to) be the point where this maximum 
is realized. Since w(xo,to) > 1, the point (xo,to) must belong to the support of n. Thus |xo| < R. 
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Since w achieves its maximum at (xo,to), then at that point we have w t > and Vw = 0. 
Therefore, u t > —m! r\ — m rjt + £o an d Vu = -mVr). We will now show that the equation for u 
cannot hold at {xQ,to). 

We start by estimating ut(xa,to). 

u t {x ,t ) > -m'(t )r](x Q , t ) - m(t )d t r](xo,t ) + e ^ 
= -m'(t )r](xo,to) + m(to)^|V?7(a;o,to)| 9 + £o 

We also note that 

|Vtt(a;o,to)| = m(t)\Vr}(x ,t )\. (3.5) 

The measure of the set |{a; : u(x,s) < 0} n B\\ is used in the estimation of (— A) s u. The 
computations below are for fixed t — to, so we omit writing the time t in order the keep the 
computations cleaner. 

Since u + mrj attains its maximum at xq, u{xq) — u{xq + y) > — m (rj(xo) — rj{xQ + y)) for 
all y. Replacing this inequality in the formula for (— A) s u we can easily obtain (— A) s u(xq) > 
— m (— A) s rj(xo). However, we can improve this estimate in terms of the value \{u < 0}nSi| (this 
is the key idea from [11]). 

Let y € R n be such that u{xq +y) < 0. We estimate u(xq) + m r/(xo) — u(xq + y) —m i](xq + y) 
using that u(xo) + mrj(xo) = w(xo, to) + £o(l + to) > 1. 

u(xq) + m rj(xo) — u(xq + y) — m rj(xo + y) > — m+1 

We choose cq small so that m < 1/2 and 

u(x ) + mr](xo) - u(x + y) - mr](xo + y) > - (3.6) 
Now we estimate (— A) s u{xq, to), we start writing the integral 
f A ss f « \ f u(x Q ) -u(x + y) , 

We estimate u(xq) — u(xo, y) by below with r](xo + y) — r](xo) except at those points where x$ + y 
is in the good set G := {u < 0} fl B\ where we use (3.6) 

> -m{t )(-A) s r){x ,t ) 
(u + mr))(x ) -(u + mr))(x + y) 

|y|n+2s a V 

>-m(t )(-Ayr 1 (x ,t Q )+ [ r^dy 

J g\b t \v\ 

> -m{t )(-A) s r,{x ,t ) + c \G\ B r \ 

for some universal constant cq (this is how Co is chosen in (3.3)). 

We consider two cases and obtain a contradiction in both. Either ?y(a;o,to) > Pi or v( x Oito) < 

Pi- 

Let us start with the second. If rj(xo, to) < pi, then (— A) s 7](xq, to) < 0, then 

(-A) s u{x ,t ) >c \G\B r \ (3.7) 
Replacing (3.4), (3.5) and (3.7) into (3.2) we obtain 

e > -m'(t )r)(x ,to) +s Q + c \G\ 
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but for any C\ > this contradicts (3.3). 

Let us now analyze the case rj{xo 1 to) > Since 77 is a smooth, compactly supported function, 
there is some constant C, such that |(— A) s ?y| < C. Then we have the bound 

(-A) s u(x ,t ) > -m(to)(-A) s n{xo,t ) + c \G\ > -Cm{t ) - c \G \ B r \ 

Therefore, replacing in (3.2), we obtain 

£0 > -m'{t )ri(xo, t ) + s + c \G\ - Cm(t ) 

and we have 

-m'{t )ri(xo, t ) - Cm{t ) + c \G\ < 0. 
We replace the value of m'(t ) in the above inequality using (3.3) 

(Ci»/(xo,*o) - C)m(t ) +c (l - »7(a?o,*o))|G| < 0. 
Recalling that 77(0:0, io) > Pi, we arrive at a contradiction if C\ is chosen large enough. □ 

The following Lemma is a crucial step in the proof of Holder continuity. It says that if a solution 
of the equation in the unit cylinder Q\ = B\ x [—1,0] has oscillation one, then its oscillation in a 
smaller cylinder Q r is less than a fixed constant (1 — 9). This type of Lemmas are very common 
when proving Holder regularity results. Sometimes they are called growth lemmas, since they 
say that the oscillation in a cylinder of a solution u necessarily grows when considering a larger 
cylinder. For nonlocal equations one needs to add an extra condition giving a bound to the 
oscillation of the function u outside Q\. 

Lemma 3.3 (diminish of oscillation) . Assume that u satisfies the following equation in Q\ 

u t -b- \7u+{-A) s u = f 

Where b is a bounded vector field. There are universal constants 9 > 0, a > and Eq > 
(depending on \ \b\\L°°, the dimension n and s) such that if 

II/IIl°°(q) < £0 

oscw < 1 
Oi 

osc u < \M/r\ a for all M > 1. 

B M X[-1,0] 

where r is a constant depending on ||6||r°o and s. Then 

osc it < (1 — 9) 

Qr ~ 

Proof. Let us choose a constant R > 1 so that R = 2 + ||&||z,°° if s > 1/2 or R = R 1 ^ 23 \\b\\ LO o + 2 
if s < 1/2. 

We consider the rescaled version of u: 

u(x,t) = u(x/(2R),t/(2R) 2s ). 

The function u satisfies the equation 

u t + (2R) 1 - 2s b ■ Vu + {-A) s u = (2R)- 2s f. 



In B2R x (— 2R,0}. We want to apply Lemma 3.1 to u minus a constant. Note that either with 
s e [1/2, 1) or s £ (0,1/2) the vector field {2R) 1 ' 2s b is bounded by A = (2i?) 1 ~ 2s \\b\\ Lao . 

Let to = (supg i u + infQj u)/2 be the median value of u in Q\ (and also the median value of 
u in Q 2 r). 

The function u will stay either above to of below to in half of the points in B\ x [—2, —1] (in 
measure). More precisely, either {u < m}r\(Bi x [-2, -1]) > |Bi|/2 or {u > m}n{B 1 x[-2, -1]) > 
|-Bi|/2. Let us assume the former, otherwise we can repeat the proof with (supg 1 u — u) instead 
of u. 

We will conclude the proof as soon as we can apply Lemma 3.1 to 2(u—m). The only hypothesis 
we are missing is that 2{u — to) is not bounded above by 1 outside B2R. So we have to consider 
v = min(l, 2(u — to)) instead, and estimate the error in the right hand side of the equation. We 
prove that if a is small enough, then v satisfies 

v t - A\Vv\ + (-A) s v < e (3.8) 

for a small £0 so that we can apply Lemma 3.1. 

Note that inside Q2R, u < 1, thus v = 2{u — m). The error in the equation in Qr comes only 
from the tails of the integrals in the computation of ( — A) s v. Indeed v t = u t and Vu = Vv. Let 
us estimate (— A) s i; — (A) s u in Qr. 

We choose eo smaller than e"o (the small constant from Lemma 3.1) such that 

(-A) s v(x,t)-(-AYu(x,t)= f (^ x + y)- snpu )+^y 

j x+ y£B 2R Qi \y\ 



jx+v&bot, \y\ 



>x+y<£B 2R W 

if a is chosen small enough. 

Thus, v satisfies (3.8). We can apply Lemma 3.1 to v and conclude the proof with r := 
l/(2R). □ 

Remark 3.4. Note that the diminish of oscillation lemma does not require the vector field b to 
be Holder continuous but only bounded. The Holder continuity of b will be required to prove the 
Holder continuity of the solution u if s < 1/2 because Lemma 3.3 needs to be a applied at all 
scales. Thus, the regularity of b compensates the bad scaling of the equation. This observation 
helps obtain some regularity results for some nonlinear supercitical equations in [10] and [5]. 



4 Holder regularity 

In this section we prove our main result about Holder continuity. Note that the result stated in 
this section requires the function to solve the equation only in a cylinder in order to have Holder 
continuity in a smaller cylinder. This is a local result, from which the global result of Theorem 
1.1 immediately follows. 

For the case s G [1/2, 1), the proof is a simple adaptation of the proof in [11]. For s € (0, 1/2) 
we need to make a change of variables at the beginning following the flow of vector field b in order 
to obtain the extra cancellation which allows us to use the hypothesis b € C 1_2s . 

Theorem 4.1. Let u be a bounded function in R" x [—1,0] which satisfies the following equation 
m Bi x [-1,0] 

u t ~b-Vu+{-A) s u = f (4.1) 
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There there is an a > such that the function u is C a at (0, 0). Moreover we have the estimate 
|u(a?,*)-u(0,0)| < C(|a;| Q + |i| Q/(2s) )(|M| ioc + ||/|| ioc ) for every x e K" and -1 < t < (4.2) 
TTie constants C and a depend on s, and the dimension n only. 

Proof. For any (xo,to), we consider the normalized function 

vM = nu\\ L J\\f\\ L ,/e< x ^ 

where eo is the constant from Lemma 3.3. 

We prove the C a estimate (4.2) byprovingaC" estimate for v at (0,0). Note that oscR» x r_i.oi v fs 
1. Moreover, v is also a solution of an equation like (4.1) and the corresponding right hand side / 
has L°° norm less than eo 

In the supercitical case s < 1/2, we need to compensate the bad scaling of the drift term with 
the required Holder regularity of b. In order to take advantage of this, it is convenient to have 
6(0, t) = 0, so that the Holder continuity translates into a decay in b close to the origin. Let us 
make a change of variables v(x, t) — v(x + A(t),t), where A(t) is a solution of the ODE: 

A(0) = 

A'(t) = b(A(t),0). 

Note that the assumption that b is Holder continuous does not guarantee the uniqueness of the 
curve A(t) following the flow of b. However, just from the continuity of b there exists one solution 
A. The uniqueness of this solution is not relevant for the proof below. 
Thus, v satisfies the equation 

v t + b ■ Vv + (- A) s v = f(x + A(t),t). 

Where b(x,t) = b(x + A(t),t) — b(A(t),t) is a Holder continuous vector field such that \b(x,t)\ < 
C\x\ 1 ~ 2s . Note that (x, t) i-> (a + A(t),t) is a Lipschitz change of variables, so the Holder regularity 
of v and v are the same. From this point on we will write v and b but we mean v and b if we 
are considering the supercritical case. For the critical or subcritical case, no change of variables is 
necessary. 

Let r € (0, 1) be as in Lemma 3.3. We will prove the following exponential decay of the 
oscillation in cylinders for some a > to be chosen later. 

oscw<r Qfe for k = 0,1,2,3,... (4.3) 
Q rk ~ 

The Holder estimate follows immediately from (4.3). 

We prove (4.3) by induction in k. Indeed, for k = it is known to be true from the hypothesis 
osc 8"x[-i,o] v ^ 1- Let us assume that it holds until certain value of k. We scale again by 
considering 

w(x,t) = r- aak v(r k x,r 2sk t). 
As it was discussed in section 2, w satisfies the equation 

w t - r ka{2s - 1] b(r k x, r 2sk t) ■ Vw + (-A) s w = r k{2s ~ a) f 

Since we are choosing a < 2s, \\ r k ( 2s - a ) f\\ Loo < e . 

In the case 2s > 1, r kal , 2s ~ 1 ) < \^ so the L°° norm of the drift term stays bounded. In the 
case 2s < 1, we use the Holder decay of b to get a uniform L°° bound for r feQ ( 2s - 1 ) b(r k x, r 2sk t) 
for \x\ < 1. 
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From the inductive hypothesis, we also have 



oscw < 1 (4.4) 
Qi 

osc to < r~ a ° k for k = 1,2, ... therefore osc w(x,t) < \M/r\ a ° for any M > 1. (4.5) 

Q r - fc " B M X[-1,0] V ^ 

If a is small enough, we can apply Lemma 3.3 to obtain oscq,, to < 1 — 9. So, if a is chosen 
smaller than the one of Lemma (3.3) and also so that 1 — 6 < r a , we have oscQ r w < r a , which 
means oscq fc+1 v < r Q ( fe + 1 ) ; which finished the proof by induction. 

We have finished the proof of (4.3) for some small a > depending on s, ||&||l°° and n. □ 



5 Weak solutions and generalizations 

In order to keep this article short and simple, we did not state the theorems in their most general 
form. We analyze in this last section the possible generalizations and the applicability of the result 
to weak solutions of (1.1). 



5.1 Viscosity solutions 

When b is a continuous vector field, the solution u of (1.1) can be understood in the viscosity 
sense. Note that in the supercitical case s < 1/2, we assume that b is Holder continuous. 

By definition, a lower semicontinuous function it is a supersolution of (1.1) in a set if every 
time there is a smooth function ip such that for some (xa,to) £ Q and p > 0, ip(xo,to) = u(xo,to) 
and tp(x,t) < u(x,t) for all points {x,t) such that \x — xq\ < p and to — p < t < to (we say <p 
touches u from below at (xo,to)), then the function 

I <p(x, t) if \x — xq\ < p and \t — to\ < p 
v(x,t) = < 

lf(a;,i) otherwise 

satisfies the inequality v t (xo,to) +b(xo,to) ■ Vv(xo,to) + (—A) s v(xo,to) > f(xo,to). 

Correspondingly, an upper semicontinuous function u is a subsolution if every time if touches 
u from above, the above constructed function v satisfies the opposite inequality. A continuous 
function u is said to be a solution when it is a subsolution and a supersolution at the same time. For 
more details on the definition and properties of viscosity solutions for intcgro-diffcrcntial equations 
see [1] or [3]. 

The definition of viscosity solution is appropriate for elliptic or parabolic equations that satisfy 
a comparison principle. As we can see in the definition, it requires that the operator can be 
evaluated in smooth test functions. Hence, the vector field b needs to be continuous. On the other 
hand, the definition of viscosity solutions does not require the equation to be linear. Thus, in the 
critical and subcritical case, the equation (1.1) can be replaced by the couple of inequalities 

u t +A\Vu\ + (-A) s it > B 
u t -A\Vu\ + (-A) s u < B 

where A = ||6|| ioo and B = \\f\\ Lao . 

The proofs presented in this article follow naturally the viscosity solution spirit. Indeed, the 
adaptation of the proof of Lemma 3.1 is straight forward since the function w constructed in that 
proof acts as the smooth test function touching u from above. 
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5.2 The vanishing viscosity method 

The results of this paper can be used to study the regularity of solutions to nonlinear equations. 
In [11], the Holder estimates were applied to the derivative of the solution to a Hamilton- Jacobi 
equation with fractional diffusion. The author is also planning to apply the main theorem in this 
paper to study the regularity of conservation laws with fractional diffusion. For these nonlinear 
problems, one can use the vanishing viscosity approximation which makes the solution smooth 
(depending on the artificial viscosity) and then obtain some estimates that are independent of the 
viscosity coefficient. In this subsection, we show that our results apply in the vanishing viscosity 
setting. 

For e > 0, we can approximate the solution to (1.1) with the solution to 

u\ + b ■ Vu e + (-A) s u E - sAu £ = f. (5.1) 

Following the arguments in this article, one can prove that the solution u £ is Holder continuous 
with an a priori estimate independent of e. 

The proof of Lemma 3.1 applies to (5.1) as long as e is bounded above by some constant C (for 
example by one). Let us restate the lemma here so that it is suitable for the vanishing viscosity 
method. 

Lemma 5.1 (point estimate). Let u be a C 2 function, u < 1 in W 1 x [—2,0] which satisfies the 
following inequality in Br x [—2,0]. 

u t - A\Vu\ + (-A) s u-eAu<e (5.2) 

Where R — 2 + A and e < 1. Assume also that 

|{u<0}n(Bi x [-2,-l])| >fx. 

Then, if so is small enough there is a 9 > such that u < 1 — 9 in B\ x [—1,0]. 
( the maximal value of So as well as the value of 9 depend only on s and n ) 

The proof is the same as the proof of Lemma 3.1 except that at the moment when we find 
a contradiction for the equation at the maximum point of w we also need to estimate eAu. We 
would have — sAu{xo,to) > em(t)A x rj{xQ,to), and then we would argue that if rj is small enough 
An > or if r\ is large, the term Cm(t)r] controls em(t)A x r](xo, to) because e < 1. 

Lemma 3.1 is used at every scale in the iterative improvement of oscillation to prove Theorem 
4.1. When scaling u to large values of k in that proof the coefficient in front of the Laplacian 
may grow to values larger than one. The result of Theorem 4.1 is still valid for values of e larger 
than zero. The proof would follow the following lines. First we do the same iteration to improve 
the oscillation in cylinders using Lemma 5.1. At some scale k, the coefficient e in front of the 
Laplacian becomes larger than one for the first time. At that point we cannot apply Lemma 5.1 
anymore, but the equation is just a lower order perturbation of the heat equation, so it has Holder 
estimates and the decay of the oscillation in the following cylinders follows because of that. 

Alternatively, we can apply Lemma 5.1 iteratively until we reach some scale k that will be finer 
the smaller e is. So, for every e > 0, we can prove that \u(x) — u(y)\ < C\x — y\ a for \x — y\ > r~ ka 
just by iterating Lemma 5.1 as in the proof of Theorem 4.1. When we take e — > 0, this estimate 
will hold for all x and y since ko —> oo as e — > 0. This is the approach that was used in [5]. 

5.3 Integro-differential diffusions 

The fractional Laplacian (— A) s in (1.1) could be replaced by general integro-differential operators 
of order 2s. These operators have the form 

Iu(x)= f (u{x)-u{x + y))^0- s Ay 

jr" \v\ 
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The kernel k(x, y) is assumed to be symmetric: k(x, y) = k(x, —y) and bounded below and above: 
A < k{x,y) < A. No modulus of continuity with respect to x is necessary. 

We now explain how to understand weak solutions of an equation of the form 



lit + b ■ Vu + 




(5.3) 



for some kernel k, some bounded right hand side / and some vector field b (either bounded or 
C l ~ 2s depending on s). If s 6 (0,1/2), the vector field b is continuous and the equation (5.3) 
implies that the following two inequalities hold: 



These operators were first defined in [9]. Even though it seems hard to make sense of (5.3) in any 
weak sense, the two inequalities above make perfect sense in the viscosity sense, and our Holder 
estimates depend on those two inequalities only. 

In the case s € [1/2, 1) the vector field b is discontinuous, so the inequalities above do not make 
sense in the viscosity sense. However, we can replace the term b ■ Vit by ±A|Vu| respectively, 
where A = \\b\\ Lao , in which case the viscosity solution definition can be applied. 

In [11], the Holder estimates were proved for a function u satisfying both inequalities above in 
the viscosity sense in the case s = 1/2. 
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